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Recently, Huang, Wu and Florin posted a Comment [l| 
on our preprint 2] describing nonequilibrium circulation 
of a colloidal sphere trapped in a optical tweezer. The 
Comment suggests that evidence for toroidal probabil- 
ity currents obtained from experiments and simulations 
in should be considered inconclusive. The authors' 
concerns are based on two claims: (1) that Brownian dy- 
namics simulations of the trapped particle's motions re- 
veal no statistically significant circulation, and (2) that 
a realistic description of the radiation pressure acting on 
the trapped sphere is inconsistent with the motion de- 
scribed in Ref. . In this Reply, we demonstrate both of 
these claims to be incorrect, and thus the original results 
and conclusions in Ref. to be still valid. 

The system, shown schematically in Fig. [l] consists 
of a single colloidal sphere trapped in a conventional 
optical tweezer formed by bringing a beam of light to 
a diffraction-limited focus 0]. In Ref. [2|, we modeled 
the trap as a radially symmetric harmonic well within 
which radiation pressure exerts an additional force di- 
rected along z: 

Fo(r) = -fcr + /iexp(^-^^ 2. (1) 

The particle's position r is measured from the trap's fo- 
cus, k is the trap's stiffness, /i sets the scale for the ra- 
diation pressure, and a is the effective range over which 
the focused light exerts forces on the particle. We assume 
that the particle is stably trapped, so that e = fi/{ka) 
may be treated as a small parameter. 

Were these the only forces acting on the sphere, the 
particle would come to a stable mechanical equilibrium 
at a distance zq ~ ecr downstream of the focus. The par- 
ticle also is acted on by random thermal forces, however, 
which displace it away from its equilibrium point. Refer- 
ence Q demonstrates analytically that the second term 
in Eq. ([T]) biases the trapped particle's thermal fluctua- 
tions in favor of toroidal circulation in the sense depicted 
in Fig. m 

Such a bias toward nonequilibrium circulation would 
occur in any model for the radiation pressure whose curl 
does not vanish. The particular choice in Eq. ([T]) facili- 
tates an analytic treatment of the effect. On this basis, 
we have claimed Q that a particle trapped in an op- 
tical tweezer does come to equilibrium, but rather acts 
as a Brownian motor 0j S @] j with the nonconservative 
component of F[){r) biasing thermal fluctuations in the 
manner of a thermal ratchet. 

The authors of Ref. [l| do not call this result into ques- 
tion, but rather claim that it is not conclusively demon- 
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FIG. 1: The probability P{0) compiled from 1000 Brownian 
dynamics simulations for a particle to circulate in an optical 
tweezer at an average rate Q. over the course of 1000 s. Inset: 
the experimental geometry, with a colloidal sphere localized 
near the focus of a beam of light propagating in the z direc- 
tion. The broad arrow indicates circulation in the positive 
direction. 

strated by the simulations and experiments presented in 
Ref. 0. To observe the predicted circulatory bias in 
a trajectory r{t) discretely sampled over time intervals 
T — 1/30 s, we introduced a measure of the mean circu- 
lation rate 

27T ^{{p-{p)Y){{z-{z)Yy 

where r = [p, cp, z) is measured in cylindrical coordinates 
centered on the trap's focal point, with z pointing along 
the optical axis. The predicted nonequilibrium circula- 
tion corresponds to clockwise rotation in the (p, z) plane, 
and to positive values of ri(t). 

Huang et al. claim [l| that Brownian dynamics simu- 
lations corresponding to the experimental conditions in 
Ref. show no statistically significant trend in fl{t), 
and thus no evidence for circulation. Our numerical sim- 
ulations, whose results are presented in Fig. [1] demon- 
strate this claim to be incorrect. Here, we have performed 
fourth-order Runge-Kutta integration of a particle's tra- 
jectory evolving according to the Langevin equation 

7f(i) = Fo(r) + Fi(<), (3) 

where 7 = Girrja is the Stokes drag coefficient for a sphere 
of radius a moving through a fiuid of viscosity 77, and 
Fi{t) is a zero-mean stochastic force whose variance is 
the thermal energy scale. Although the simulations in 
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[2| were performed with the isotropic model in Eq. ((T)), 
the authors of Ref. Q generahze the harmonic restoring 
force to account for different trap stiffnesses in the three 
Cartesian directions, using = 0.467 pN//im, ky — 
0.4 pN//^m and kz — 0.08 pN/^m. In responding to 
their criticism, we adopt the same anisotropic force law in 
simulating the motions of a sphere of radius a = 1.1 /im in 
a trap of width a = a. Following Refs. [H and 0], we also 
set e — 0.1 and adopted time steps of 10"'* s. Taking the 
suspending medium to be water at room temperature, 
ry = 10-3 Pas. 

Figure [T] shows the distribution of mean circulation 
rates, = {fl{t)), obtained from 1000 independent runs, 
each of 1000 s duration. As Huang et al. point out [ij, 
individual realizations can display either positive or neg- 
ative circulation. Contrary to their assertion, however, 
these variations do not occur with equal probability. Out 
of 1000 realizations, only 131 showed negative circulation. 
A similarly small proportion of retrograde circulation is 
observed experimentally. The ensemble-averaged circu- 
lation rate fl = 0.08 Hz agrees quantitatively both with 
the experimental results and also with the analytic pre- 
dictions presented in Ref. 0. 

Although the principal claim by Huang et al. is thus 
shown to be incorrect, their Comment raises the valid 
point that the scattering force experienced by a colloidal 
sphere in a real optical trap is likely to be more compli- 
cated than the idealized model in Eq. ([T]), particularly 
for spheres larger than the wavelength of light. The de- 
tailed form of Fq is less important, however, than the 
presence of a rotational component, V x Fq ^ 0, for bi- 
asing the system out of equilibrium. It is this rotational 
component that breaks the spatiotemporal symmetry of 
the particle's fluctuations to create a net flux in its prob- 
ability density 6]. The particular form in Eq. ([T]) was 
selected more for its analytic tractability than for its ac- 
curacy as a model for radiation pressure in optical traps. 

The form for the scattering force presented in Ref. [Ij 
also has a rotational component, and so will give rise to 
circulation in the particle's trajectory. Unlike the scat- 
tering force in Eq. |T]), which is peaked on the optical 
axis, the ray-optics calculation in Ref. [H increases with 
distance from the optical axis, and so would induce ret- 
rograde circulation. This observation raises the interest- 
ing point that optically trapped particles' behavior may 
be more complicated than is predicted by the idealized 
model in Eq. ([T]). 

If this model for the radiation pressure were relevant 
to the experiments in Ref. i^, then the measured trajec- 



tories also should have displayed retrograde circulation. 
Huang et al. suggest that the discrepancy can be ascribed 
to insufflcient statistics in the experimental analysis. We 
argue instead that the result for the scattering force pre- 
sented in Fig. 2 of Ref. [H reflects only the z component 
of the optical force that a particle would experience at 
the optical tweezer's focal point. In fact, the silica sphere 
is more than twice as dense as the water in which it is 
suspended, and so settles roughly 2 ^m below the focal 
point. In this region of the beam, the total optical force 
computed by a fully vectorial theory [3] has a uniformly 
positive curl. Consequently, the particle should undergo 
positive circulation, as reported. This is not to say that 
Eq. ^ is an accurate representation for the optical forces 
experienced by the sphere, but rather that the form pro- 
posed in Fig. 2 of Ref. [l[ is not. 

In conclusion, we have demonstrated that the concerns 
raised by Huang, Wu and Florin in Ref. 1] can be as- 
cribed to inadequate statistical analysis of their simu- 
lations and to an incomplete analysis of the scattering 
force acting on optically trapped spheres. The results 
and conclusions presented in Ref. [2] therefore remain 
unchanged. 
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